Quantum Computation as Geometry 
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^ ■ Quantum computers hold great promise, but it remains a challenge to find 

cn : 

^ ; efficient quantum circuits that solve interesting computational problems. We 

o ■ 

; show that finding optimal quantum circuits is essentially equivalent to finding 

^1 the shortest path between two points in a certain curved geometry. By recast- 

§ ■ ing the problem of finding quantum circuits as a geometric problem, we open 

^ : 

up the possibility of using the mathematical techniques of Riemannian geom- 

> ' 

■ etry to suggest new quantum algorithms, or to prove limitations on the power 

^' 

I of quantum computers. 

Quantum computers have the potential to efficiently solve problems considered intractable 
on conventional classical computers, the most famous example being Shor's algorithm (jij) for 
finding the prime factors of an integer. Despite this great promise, as yet there is no general 
method for constructing good quantum algorithms, and very little is known about the potential 
power (or limitations) of quantum computers. 

A quantum computation is usually described as a sequence of logical gates, each coupling 
only a small number of qubits. The sequence of gates determines a unitary evolution U per- 
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formed by the computer. The difficuky of performing the computation is characterized by the 
number of gates used by the algorithm, which is said to be efficient if the number of gates re- 
quired grows only polynomially with the size of problem (e.g. with the number of digits in the 
number to be factored in the case of Shor's factoring algorithm). 

We develop an alternate approach to understanding the difficulty of implementing a uni- 
tary operation U . We suppose that U is generated by some time-dependent Hamiltonian H{t) 
according to the Schrodinger equation dU /dt = —iHU, with the requirement that at an appro- 
priate final time U{tf) = U. We characterize the difficulty of the computation by imposing 
a cost F{H{t)) on the Hamiltonian control, H{t). Following we choose a cost function 
on H{t) that defines a Riemannian geometry on the space of unitary operations. Finding the 
optimal control function H{t) for synthesizing a desired unitary U then corresponds to finding 
minimal geodesies of the Riemannian geometry. 

We will show that the minimal geodesic distance between the identity operation and U is 
essentially equivalent to the number of gates required to synthesize U. This result extends 
the work in (2 ), where it was shown that the minimal distance provides a lower bound on the 
number of gates required to synthesize U. 

The power and interest of our result is that it allows the tools of Riemannian geometry to 
be applied to understand quantum computation. In particular, we can use a powerful tool — 
the calculus of variations — to find the geodesies of the space. Just as in general relativity, this 
calculus can be used to derive the geodesic equation, a "force law" whereby the local shape of 
space tells us how to move in order to follow the geodesies of the manifold. 

Intuitively, our results show that the optimal way of solving any computational problem is to 
"fall freely" along the minimal geodesic curve connecting the identity operation to the desired 
operation, with the motion determined entirely by the local "shape" of the space. To appreciate 
how striking this is, consider that once an initial position and velocity are set, the remainder of 
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the geodesic is completely determined by the geodesic equation. This is in contrast with the 
usual case in circuit design, either classical or quantum, where being given part of an optimal 
circuit does not obviously assist in the design of the rest of the circuit. Geodesic analysis thus 
offers a potentially powerful approach to the analysis of quantum computation. However, a 
caveat to this optimism is that although we know the initial position is the identity operation, 
we still need to determine the initial velocity in order to find the minimal geodesic, and this is 
not in general an easy problem. 

Our results can also be viewed as showing that the problem of finding minimal quantum 
circuits is equivalent to a problem in geometric control theory ©, which has had great success 
in using techniques from the calculus of variations and Riemannian geometry to solve optimal 
control problems. For example, Khaneja et al (4j (c.f. also have used geometric tech- 

niques to analyse the minimal time cost of synthesizing two-qubit unitary operations using a 
fixed two-qubit control Hamiltonian, and fast local control. 

In order to choose a cost function on the control Hamiltonian H(t) we first write H{t) in 
terms of the Pauli operator expansion H = J2'^ h^a + X)" h^jU, where: (1) in the first sum a 
ranges over all possible one- and two-body interactions, that is all products of either one or two 
Pauli matrices acting on n qubits; (2) in the second sum a ranges over all other tensor products 
of Pauli matrices and the identity; and (3) the h„ are real coefficients. We then define a measure 
of the cost of applying a particular Hamiltonian during synthesis of a desired unitary operation 



F{H) 
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The parameter p is a penalty paid for applying three- and more-body terms; later we will choose 
p to be large, in order to suppress such terms Q. 

This definition of control cost leads us to a natural notion of distance in SU (2") . A curve \U] 
between the identity operation / and the desired operation f/ is a smooth function U : [0, t/] — > 



SU{2"-) such that U{0) = I and U{tf) = U. The length of this curve can then be defined by 
the total cost of synthesizing the Hamiltonian that generates evolution along the curve: 



Since (i([f/]) is invariant with respect to different parameterizations of [U] we can always 
rescale the Hamiltonian H{t) such that F(H(t)) = 1 and the desired unitary U is generated at 
time tf = d{[U]). From now on we assume that we are working with such normalized curves. 
Finally, the distance d{I, U) between / and U is defined to be the minimum of rf([f/]) over all 
curves \U] connecting / and U . 

We will show that for any family of unitaries U (implicitly, U is indexed by the number of 
qubits, n) there is a quantum circuit containing a number of gates polynomial in U) that ap- 
proximates U to high accuracy. In other words, if the distance d{I, U) scales polynomially with 
n for some family of unitary operations, then it is possible to find a polynomial-size quantum 
circuit for that family of unitary operations. Conversely, the metric we construct also has the 
property, proved in (0), that up to a constant factor the distance d{I, U) is a lower bound on the 
number of one- and two-qubit quantum gates required to exactly synthesize U . Consequently, 
the distance d{I, U) is a good measure of the difficulty of implementing the operation [/ on a 
quantum computer. 

The function F{H) specified by Eq.[l]can be thought of as the norm associated to a (right 
invariant) Riemannian metric whose metric tensor g has components: 



These components are written with respect to a basis for the local tangent space corresponding 
to the Pauli expansion coefficients h„. The distance d{I, U) is equal to the minimal length solu- 
tion to the geodesic equation, which may be written (9) as {dH/dt, K) = i{H, [H, K]). In this 
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' if a ^ r 
1 if (J = r and a is one- or two-body 
if cr = r and a is three- or more-body. 



(3) 



expression, (-, ■) is the inner product on the tangent space sm(2") defined by the metric compo- 
nents of Eq.|3l and K is an arbitrary operator. For our particular choice of metric components, 
this geodesic equation may be rewritten as: 

pIK = i^plKh[a,T], (4) 

r 

where hy„^r\ = tr(iJ [cr, r])/2". A particular class of solutions to this equation was studied 
in (|2l), but understanding the general behaviour of the geodesies remains a problem for future 
research (I7DI) . We note that there are powerful tools in Riemannian geometry (see, e.g., (Iiilli2l) ) 
available for the study of minimal length geodesies. 

Our goal is to use the optimal control Hamiltonian H{t) to explicitly construct a quantum 
circuit containing a number of gates polynomial in d{I, U), and which approximates U closely. 
The construction combines three main ideas, which we express through three separate lemmas, 
before combining them to obtain the result (Fig. 1). 

The first lemma shows that the error that arises by simply ignoring the many-body inter- 
actions in H(t) can be made small by choosing the penalty p appropriately. We define Hp to 
be the projected Hamiltonian formed by deleting all three- and more-body terms in the Pauli 
expansion. Then the following result is proved in the supporting online materials. 

Lemma 1: Let Hp{t) be the projected Hamiltonian obtained from a Hamiltonian H(t) 
generating a unitary U . Let f/p be the corresponding unitary generated by Hp(t). Then 

||f;_t,HI<«. (5) 

P 

where || ■ || is the operator norm (TT^ . and p is the penalty parameter appearing in the defi- 
nition of the metric. Thus, by choosing p sufficiently large, say p = 4", we can ensure that 
\\U - Up\\ < rf([f/])/2". 

Motivated by the preceding lemma, we change our aim from accurately synthesizing U to 
accurately synthesizing Up. To do this, we break the evolution according to Hp(t) up into 

5 



Figure 1 : Schematic of the three steps used to construct a quantum circuit approximating the 
unitary operation U . The circuit is of size polynomial in the distance d{I, U) between the 
identity and U . First we project the Hamiltonian H{t) for the minimal geodesic path onto one- 
and two-qubit terms, giving Hp(t). By choosing the penalty p large enough (p = 4") we ensure 
the error in this approximation is small, ei < d{I, [/)/2". Next we break up the evolution 
according to Hp{t) into small time steps of size A = U)/N, and approximate with a 
constant mean Hamiltonian Hp over each step. Finally we approximate evolution according to 
the constant mean Hamiltonian over each step by a sequence of one- and two-qubit quantum 
gates. The total errors, 62 and 63, introduced by these approximations can be made smaller than 
any desired constant by choosing the step size A sufficiently small, A = 0{1/ {n^d{I , U))). 
In total we need 0{n^d{I, U)^) quantum gates to approximate U to within some constant error 
which can be made arbitrarily small. 
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many small intervals, each of length A. The next lemma shows that evolution according to the 
time-dependent Hamiltonian Hp{t) over such a small time interval can always be accurately 
simulated by a constant mean Hamiltonian, which we denote Hp. 

Lemma 2: Let U be an n-qubit unitary generated by applying a time-dependent Hamilto- 
nian H{t) satisfying \\H(t) \\ < c over a time interval [0, A]. Then defining the mean Hamilto- 
nian H = ^ dt H{t) we have: 

\\U - exp(-ii7A)|| < 2(e"^ - 1 - cA) = ©(c^A^). (6) 

The proof of this lemma is based on the Dyson operator expansion and is presented in 
the appendix. To apply this lemma to Hp{t), note that elementary norm inequalities and the 
observation F(i7p(t)) < 1 imply that (dl \\Hp{t)\\ < ^nF{Hp{t)) < -^n. Lemma 2 
implies that over a time interval A we have: 



2A2^ 



||[/^ -exp(-2i?^A)|| < 2 l^e^/^"^ - + -^n/^ = Oiji^A 

where Up is the evolution generated by Hp(t) over the time interval A, and Hp is the corre- 
sponding mean Hamiltonian. 

Our third and final lemma shows that evolution according to a time-independent Hamil- 
tonian H containing only one- and two-body terms can be very accurately simulated using a 
number of quantum gates that is not too large. 

Lemma 3: Suppose H is an n-qubit two-body Hamiltonian whose Pauli expansion coeffi- 
cients satisfy \h(j\ < 1. Then there exists a unitary Ua, satisfying 

||e-*^^-f/^|| < Can^A^, (8) 

that can be synthesized using at most c\r? j l\ one- and two-qubit gates, where c\ and C2 are 
constants. 
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This result follows from standard procedures for simulating quantum evolutions using quan- 
tum gates (see, e.g., Chapter 4 of (175b ). and is proved in the appendix. Note that the average 
Hamiltonian Hp provided by Lemma 2 satisfies the assumptions of Lemma 3, since the Pauli 
expansion coefficients of Hp{t) satisfy \h^\ < 1 for all times. 

To integrate Lemmas 1-3, suppose H{t) is the time-dependent normalized Hamiltonian gen- 
erating the minimal geodesic of length d{I,U). Let Hp(t) be the corresponding projected 
Hamiltonian, which generates Up and satisfies ||f/ — f/p|| < d{I,U)/2'", as guaranteed by 
Lemma 1, and where we have chosen p = 4" as the penalty. Now divide the time interval 
[0, d{I, U)] up into a large number of time intervals each of length A = d(I, U)/N. Let Up 
be the unitary operation generated by Hp{t) over the jth time interval. Let f/{^ be the unitary 
operation generated by the corresponding mean Hamiltonian. Then Lemma 2 implies that: 

\\U'p - UUW < 2(e3/v^"^ - (1 + -^nA)). (9) 

Lemma 3 implies that we can synthesize a unitary operation using at most c\r? j IS. one- and 
two-qubit gates, and satisfying \ U\^ — ^All ^ C2n'^A^. 

Putting all these results together and applying the triangle inequality repeatedly, we obtain: 

\\V-Ua\\ < \\U-Up\\ + \\Up-Ua\\ (10) 

< '-^^hu^^U'J (H) 

< ^^ + j:{\m-UiA\ + \\Uh-Ui\\) (12) 

i=i 

< + 2^ (e(^/^)"^ - (l + -|,a) ) + c,d(/, U)n^ A\(13) 

Provided we choose A to scale at most as l/{n^d{I, U)), we can ensure that the error in our 
approximation Ua to U is small, while the number of gates scales as n^d{I, U)^. Summing up, 
we have the following theorem (173b : 
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Theorem: Using 0{n^d{I ,Uy) one- and two-qubit gates it is possible to synthesize a 
unitary Ua satisfying \\U — Ua\\ < c, where c is any constant, say c = 1/10. 

Our results demonstrate that, up to polynomial factors, the optimal way of generating a 
unitary operation is to move along the minimal geodesic curve connecting / and U. Since 
the length of such geodesies also provides a lower bound on the minimal number of quantum 
gates required to generate U, as shown in ©, the geometric formulation offers an alternative 
approach which may suggest efficient quantum algorithms, or provide a way of proving that a 
given algorithm is indeed optimal. 

It would, of course, be highly desirable to completely classify the geodesies of the metric 
we construct. An infinite class of such geodesies has been constructed in and shown to have 
an intriguing connection to the problem of finding the closest vector in a lattice. In future, a 
more complete classification of the geodesies could provide significant insight on the potential 
power of quantum computation. 



Proof of Lemma 1: We require three facts about the (unitarily-invariant) operator norm || ■ || 
and the cost function F{H): 

1. Suppose time-dependent hamiltonians H(t) and J(t) generate unitaries U and V , respec- 
tively, according to the time-dependent Shrodinger equation. By repeatedly applying the 
triangle inequality and the unitary invariance of || ■ ||, we obtain the inequality: 




2. If H contains only three- and more-body terms we have F{H) = p\\H\\2, where || ■ ||2 is 
the Euclidean norm with respect to the Pauli expansion coefficients. 



Appendix: 





(14) 
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3. For any H 



\H\ 



a 



<El^-l<2l|i^||2, 



where the final inequality follows by an application of the Cauchy-Schwarz inequality. 



Combining these observations we have 



d{[U]) = J dtF{H{t)) 

> J dtF{H{t) - Hp{t)) 

> j dtp\\H{t)-Hp{t)h 

> l-j dt\\H(t)-Hp{t)\\ 



> ^\\U-Up\\, 

from which the result follows. 

Proof of Lemma 2: Recall the Dyson series (see, e.g., (17^ . p. 325-326): 



(15) 
(16) 
(17) 
(18) 
(19) 



m=0 



dti / dt2 



dtmH{U)H{t2)...H{t^). 



(20) 



Note that in our finite-dimensional setting this series is always convergent. By writing out 
the power series for exp(— zffA) and canceling the 0(A°) and O(A^) terms and applying the 
triangle inequality we have 

-ii/A)™!! 



-u\\ < E 



m=2 

j-A j-h 

+ dti dt2 
Jo Jo 

< 2E 



tm — 1 



m=2 



dt^\\H{ti)H{t2)---H{t, 
2 (e^^ - 1 - cA) , 



(21) 
(22) 



where for the second line we have used the standard norm inequality ||Xy|| < ||X||||F||, the 
condition \\H{t)\\ < c and the fact that dti !o dt2--- /o"~' dU, = A^/m!. This is the 
required result. QED 
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Proof of Lemma 3: We use standard quantum simulation techniques, see e.g. Section 4.7 
of (175b . Divide the interval [0, A] up into N = Ij/S. steps of size A^. Define 

U^2 = e-^^^<^^^\-ih2^2^^ _ _ _ g-iftL<TiA2 ^23) 

where L = 0{n^) is the number of terms (all one- and two-body) in H . Each factor on the right 
hand side of this equation is an allowed one- or two-qubit quantum gate. Using, e.g., Eq. (4. 103) 
on page 208 of ( 15 ) it is straightforward to show that: 

^7^2 =e-^^^' + 0(L2A^), (24) 

and therefore through repeated applications of the triangle inequality and the unitary invariance 
of the operator norm we obtain: 

||g-iHA _ II < caiVn^A^ = csn^A^, (25) 

where C2 is a constant. Eq. (l25t shows how to approximate e~*^^ using at most Ci^^/A quan- 
tum gates, where ci is another constant, which is the desired result. QED 
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